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Statistical fluctuation problems are faced by all quantum key distribution (QKD) protocols under 
finite-key condition. Most of the current statistical fluctuation analysis methods work based on 
independent random samples, however, the precondition cant be always satisfied on account of dif¬ 
ferent choice of samples and actual parameters. As a result, proper statistical fluctuation methods 
are required to figure out this problem. Taking the after-pulse contributions into consideration, 
we give the expression of secure key rate and the mathematical model for statistical fluctuations, 
focusing on a decoy-state QKD protocol (Sci Rep. 3, 2453, 2013) with biased basis choice. On 
this basis, a classified analysis of statistical fluctuation is represented according to the mutual re¬ 
lationship between random samples. First for independent identical relations, we make a deviation 
comparison between law of large numbers and standard error analysis. Secondly, we give a sufficient 
condition that Chernoff bound achieves a better result than Hoeffdings inequality based on indepen¬ 
dent relations only. Thirdly, by constructing the proper martingale, for the first time we represent 
a stringent way to deal with statistical fluctuation issues upon dependent ones through making use 
of Azumas inequality. In numerical optimization, we show the impact oil secure key rate, the ones 
and respective deviations under various kinds of statistical fluctuation analyses. 

PACS numbers: 03.67.Dd, 42.81.Gs, 03.67.Hk 

I. INTRODUCTION 


Quantum key distribution (QKD) protocols can 
provide unconditional security for telecommunications 
based on the physical law of quantum mechanics [I-8]. 
However, in real-life QKD systems, there is a large gap 
between the ideal conditions and the practical ones. In 
order to figure out the problem, Gottesman et al. [9j 
gave an intensive study aiming at the real system im¬ 
perfections, and provided a key rate computation for¬ 
mula called GLLP. Based on GLLP formula, decoy-state 
protocols floUTl were proposed so as to improve se¬ 
cure key rate and detect photon-number-splitting (PNS) 
attack[lj, [jj|. For the first 20 years since the proposal 
of BB84 protocol [lj , QKD postprocessing procedure had 
been mainly analyzed in the asymptotic assumption that 
the length of transmission code is infinite. Nevertheless, 
in real-world situations, the two parties often commu¬ 
nicate within a limited time. As a result, many efforts 
0, 13 El-il have been spent on the finite-key effect in 
QKD postprocessing schemes. 

The security of QKD under finite-key condition was 
first considered by Mayers 0|, and soon afterwards, by 
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using standard error analysis, Ma 0 gave a detailed 
analysis concerning five kinds of statistical fluctuations. 
Based on the work in [3 H3, Wei et al. 2l| pro¬ 


posed an efficient decoy-state QKD protocol with bi¬ 
ased basis choice and adopted standard error analysis 
to tackle statistical fluctuation problems. More recently, 
Zhou et al. [3] proposed an improved statistical fluctua¬ 
tion method in measurement-device-independent (MDI) 
QKD also based on standard error analysis, which raises 
the key rate of the recent Shanghai experimental set-up 
13 by more than 600 times. 

Under universal composable (UC) [29j security which 
is accepted and considered as the most stringent secu¬ 
rity standard in classical and quantum fields, Scarani and 
Renner [3| presented the finite-key QKD security analy¬ 
sis using Law of large numbers [3, EH, HQ] . By adopting 
Law of large numbers to deal with statistical fluctua¬ 
tion issues, Cai and Scarani [18] discussed the secure key 
generation rate in their decoy-state protocol and obtain 
the shortest sending code length for generating secure 
key. Based on their research, Lucamarini et al. [241 ] used 
beta distribution and CP test [3lj to handle the finite- 
key effects and gave an experimental verification for their 
decoy-state QKD protocol. 

By taking advantage of uncertainty relationship for 
smooth entropies, Tomamichel et al. [l9| provided an up¬ 
per bound of the final key. Based on their work, through 
utilizing Hoeffding’s inequality (32| to analyze statisti- 
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cal fluctuations in decoy-state QKD protocol, Ref. [25} 
achieved a good estimation result. Since then large devi¬ 
ation theories [H], HI} have attracted lots of attentions in 
dealing with statistical fluctuation problems. Applying 
Chernoff bound @,[3^] to perform parameter estimation, 
Curty et al. demonstrated the feasibility of long-distance 
implementations of MDIQKD in finite-code length. 

In most of the above QKD postprocessing schemes 
lim2014concise [H, El [2l|, [2f| [2f| yj, El , the main 
precondition of statistical fluctuation analysis is almost 
the same, which requires an independent mutual rela¬ 
tion between random samples. Renner [36} pointed out 
that symmetry of large physical systems implies inde¬ 
pendence of subsystems, thus the independent relation¬ 
ship is reasonable as long as the systems physical prop¬ 
erties are robust under small disturbances. However, 
when it comes to after-pulse contributions, the detec¬ 
tion event might be caused by the after-pulse of the prior 
one. Many efforts have been spent on investigating after¬ 
pulse [33, HH and its effect [39t Ho} . concentrating on 
the correlations of neighbouring pulses and the security 
analysis. Somma et al. [39] demonstrated that if the in¬ 
dependent condition is broken by after-pulse contribu¬ 
tions, then the security parameters of decoy-state imple¬ 
mentations could be worse, and they also pointed out 
that the eavesdropper possesses the ability to correlate 
her PNS attack which breaks the independent condition, 
so that the corresponding security analysis would be in¬ 
valid in the above mentioned protocols. Furthermore, the 
after-pulse effect will lower the secure key rate, which is 
even worse at higher gate frequencies because of a large 
number of after-pulse-related errors j4Cj. As a result, the 
after-pulse contributions cant be neglected, and we need 
to find a stringent way to reconsider the security analy¬ 
sis, especially statistical fluctuation issues based on de¬ 
pendent random samples, which is the motivation behind 
our analysis. 

In this paper, we focus on the mathematical mod¬ 
els, classifications and deviation comparisons for various 
kinds of statistical fluctuation analyses, including stan¬ 
dard error analysis, Law of large numbers, Hoeffdings 
inequality and Chernoff bound. However, these methods 
are not the accurate ways to deal with statistical fluc¬ 
tuation issues when taking the after-pulse contribution 
into consideration. In order to solve the problem, for the 
first time we provide a new way called Azumas inequality 
[35} to deal with dependent random samples. Since the 
efficiency of a practical QKD system needs to care about 
the total quantum resource consumption, it is important 
to take the sending qubits as the total random samples. 

The article is organized as follows. Section HD gives an 
analysis of after-guise contributions based on Weis decoy- 
state protocol [2l[ . cm we set up mathematical models 
for standard error analysis and Law of large numbers 
based on independent identically distributed (i.i.d.) ran¬ 
dom samples, and present a deviation comparison. For 
independent random samples only, we give a similar anal¬ 
ysis for large deviation theories in section IIVI In section 


El we bring a stringent way to solve the statistical fluc¬ 
tuation issues on dependent samples. Simulations for se¬ 
cure key rates and respective deviations of the five kinds 
of estimation methods are given in section IVII At last, 
some conclusions are presented in section ETT1 


II. ANALYSIS OF AFTER-PULSE 
CONTRIBUTIONS 

The after-pulse effect is due to carrier traps, and the 
uncontrollable clicks will contribute to the quantum bit 
error rate(QBER). To make a detailed analysis, here we 
choose vacuum+weak decoy-state QKD protocol with bi¬ 
ased basis choice [2l|- Taking after-pulse contributions 
into consideration, we need to reconsider the expected 
detection rate, the yield of vaccum decoy-state yield and 
the observed QBER. 


A. Parameter description 


Denoting p ap as the after-pulse probability, here we 
suppose Pap = 4 x 10" 2 gl. Then the parameters men¬ 
tioned above change their expressions as follows. 


D p = Q^(l+Pa P ) = (l-e-w(l-Y 0 ))(l+Pa P ) 

D v = Q v (1 + Pap ) = (1 - e(1 - Y 0 )) (1 +Pap) 

Y 0 ' = 2p dc (l+p ap ) 
p , _ (eoV+eda-e^Hl-YoO+PapQ^A) 

“ Da 

(1) 

where Dp, D v , Yo' and E p respectively denote the over¬ 
all gain of signal pulses, the weak decoy-state pulses, 
the yield of vacuum decoy-state pulses and the observed 
QBER under after-pulse impacts; p, p dc , eo and e d rep¬ 
resent the transmittance of the channel, the background 
count rate, the background error rate and the probability 
that a photon hits the wrong detector. 

Further, we obtain the secure key rate 


R 

leakEC 

q 

Qo 


> <?{Qo + Qi[l- H (e{ z )} - leaksc} 

= fEcQpH{E p ) 

^ NqPz 
Ntotal 

= Y 0 r e~» 


( 2 ) 


where q denotes the received signal state ratio; Qo and 
Q z represent the background and the single photon gains, 
e v { is the phase error rate of single photon signal state; 
leakEC, Sec, Qp and E p separately denote the cost of 
error correction, the bilateral error correction inefficiency, 
the overall gain of signal state and quantum bit error rate 
(QBER); H (x) = —Hog 2 ( 2 ) — (1 — a,’)log 2 (l — a:)is the 
binary Shannon entropy function; N tota i , Np, Pz and Yo 
are the length of the sending qubits, the signal pulses in 
the basis, the probability that Bob chooses the basis, and 
the background count rate. 

In security analysis, we need to estimate the lower 
bound of single photon yield denoted by Yi and the up¬ 
per bound of its relevant error rate ei measured in the Z 
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basis. Their expressions are 


B. Law of large numbers 


Yl > 


and 


fiv — v 2 






ei < e 


u 

1 


E v Q v e v - e 0 Y 0 ' 
Yfv 


( 4 ) 


B. Mathematical model 

Mathematical model for the corresponding statistical 
fluctuation analysis considering after-pulse contributions 
is discussed as follows. Let X denote the observed value, 
which is obtained by detection events of m random sam¬ 
ples denoted by Xi, X 2 , ■ ■ ■ , X m , detected with the value 
1 and 0 otherwise, satisfying X = X Xj. E is the 
expected value of Y X; when m tends to infinity. 
Because of after-pulse contributions, we have the follow¬ 
ing probability relations 


Law of large numbers [13 mainly describes the stable 
equivalence relation between occurrence frequency and 
probability distribution when the test number is large 
enough. Its application precondition only requires that 
the random samples are i.i.d.. And the mathematical 
model is stated as follows. 

The only difference from the above mentioned one is 
that the observed value X is obtained by detection events 
of to i.i.d. random samples. Then for any £ > 0 and 
e > 0, a quantum version of law of large numbers yields 
the following statement 



numbers in dealing with statistical fluctuation problems 
in Wei’s protocol, we obtain the relevant deviations de¬ 
noted by £ 2 (Softs), which can be found in Appendix iBl 


Pr {X m +i — s\X m — t, X n —\ — x n —i, • • • , X\ — Xi, Xq — £ 0 } — PS — s |AT m — 7 ^ Y*v — s} , 


( 5 ) 

which shows that the random samples are Markov chain. 
Then for any and , we need to find a proper inequality 
satisfying 

Pr (|X — E\ > £) < £ (6) 

where the parameter can be seen as a function of . That 
is to say the expected value is contained in the interval 

[X-£,X + £] (7) 

with a failure probability e. 


III. ANALYSIS OF I.I.D RANDOM SAMPLES 

The purpose of this section is to make a deviation com¬ 
parison between two kinds of analysis methods based on 
i.i.d. random samples. 


A. Standard error analysis 

The original protocol adopts standard error analysis 
which is the application of central-limit theorem using 
Gaussian distribution to deal with the statistical fluctu¬ 
ation problems. The application precondition is that the 
random samples must be i.i.d. and the sum of them has 
to obey Gaussian distributions. We denote £1 (Softs) as 
the deviation when using standard error analysis, here 
Sobs refers to the respective random samples. There are 
four kinds of deviations corresponding to Wei’s protocol 
defined in Appendix [A] 


C. Deviation comparison 

Here we should point out that the quantile denoted 
by u a in standard error analysis is determined by the 
failure probability when applying the statistical fluctua¬ 
tion analysis method. Thus, before making a deviation 
comparison between the two methods, we need to set up 
a fair comparison standard first. We set a fixed failure 
probability no matter which method is used. Suppose 
that the failure probability is e = 10 10 , which is equiv¬ 
alent to taking the quantile as u a = 6.4. So we have to 
change the quantile from u a = 5 in the original protocol 
into u a = 6.4. 

Given that the length of sending qubits is long enough 
so that the sum of them satisfies Gaussian distribution. 
Now we just take Q M as an example and the respec¬ 
tive deviations are £1 (Q^) = = and £2 (Qp) = 

V NyPzQy. 

2[ln 10 lo +2 ln(+1)]~ 

N^PzQy . 2 ' 

By some reduction, it equals to compare 6.4- yjQ M and 
V46.06 + 41n(Y M p z + l). 

Since y/Q^ < 1 and 41n(A^ M p z + 1) > 0, then we can 
obtain the following inequalities 

yj 46.06 + 4 In^N^pz + 1) > v^6(06 > 6.4 > 6.4 • JQ^. 

( 9 ) 

As a result, it can be safely concluded that the rela¬ 
tionship between relevant deviations satisfies £1 < £ 2 - In 
other words, standard error analysis can achieve a better 
estimation result than law of large numbers. The simu¬ 
lation result is shown in El 
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IV. ANALYSIS OF INDEPENDENT RANDOM 
SAMPLES 

Large deviation theory [H, 33] describes the remote 
tails’ asymptotic behavior of probability distribution se¬ 
quences. That is to say, it concerns the exponential de¬ 
cline of the probability measures of certain kinds of tail 
events. It is a precise form of law of large numbers and 
the main application precondition only requires the re¬ 
lation of observed samples is independent. Here we take 
two basic theorems to deal with the statistical fluctuation 
issues concerning after-pulse contributions. 

A. Hoeffding’s inequality 

The theorem introduced first is called Hoeffding’s in¬ 
equality [32], which provides an upper bound on the prob¬ 
ability that the sum of independent random variables de¬ 
viates from its expected value. A quantum version of 
Hoeffding’s inequality yields the following statement. 

Let X denote the observed value, which is obtained by 
detection events of m independent random samples. And 
Xi are almost surely bounded, that is Pr (Xj £ [a^, bi\) = 
1, 1 < i < m.Then for any £ > 0 and e > 0, the following 
inequality holds 

Pr (| A - E [X] | > £) < = e, ( 10 ) 

here E [X] refers to the expected value of X. 

Using Hoeffding’s inequality, we get the relevant devi¬ 
ations denoted by £ 3 (Sots), which can be found in Ap¬ 
pendix 1 X 1 


Employing Cliernoff bound in Wei’s protocol, we obtain 
the corresponding deviations denoted by ^(S'ofcs) in Ap¬ 
pendix [A] 


C. Deviation comparison 

In the definition of Chernoff bound, it can be easily 
concluded that if the three conditions are not satisfied, 
then Chernoff bound changes its form into Hoeffding’s 
inequality. Here we only consider the situation that the 
three conditions in Chernoff bounds are all fulfilled, and 
the comparisons corresponding to the other results are 
shown in Appendix [B] 

For convenience, we suppose £1 = £2 = £3 = ICR 10 . 
Still we take Q M as an example and the respective devi¬ 
ations are 6 (Q M ) = yj ln R A = ln £ and 

A' = \J lnp^r, here £ 3 , A and A' refer to the de¬ 
viations when using Hoeffding’s inequality and Chernoff 
bound. 

Then if both yj-^- In if < yj^- ln ± and 

w^t ln w2 < ^ 2 w~p.- ln ^ hold ’ we can safel y sa y 

that Chernoff bound can achieves a better result than 
Hoeffding’s inequality. 

Simplify the above inequalities, we obtain the sufficient 
condition is 

Qp < 0 . 06 . ( 12 ) 


V. ANALYSIS OF DEPENDENT RANDOM 
SAMPLES 


B. Chernoff bound 


The second theorem introduced here is called Chernoff 
bound [33], also providing an upper bound on the proba¬ 
bility that the sum of independent random variables de¬ 
viates from its expected value. To apply the method in 
Wei’s protocol, we changes the form mentioned in jJI} 
into the following claim. 

Let X denote the observed value, which is obtained 
by detection events of m independent random samples, 
satisfying Pr(Xi = 1) = Pi, X = and 

a = E [X] = where E [•] denotes the mean 

value. Let x be the observed outcome of X and = 



for certain £1 > 0 . 


For certain £ 2 , 


£3 > 0 , we have that x satisfies 


x = a + <5, 


( 11 ) 


except with error probability 7 = £1 + £2 + £ 3 , where 
5 £ [—A, A']. Here £2 (£ 3 ) refers to the probability that 
x < a — A (a; > a + A'). The limitations and corre¬ 
sponding values of A and A' are given in Appendix [C] 


Concerning the worst case that the after-pulse (maybe 
caused by eavesdropper) exists in every detection event, 
the mutual relationship between adjacent samples doesnt 
satisfy the independent condition anymore. Therefore, 
we need a stringent way to deal with the problem based 
on dependent observed samples. Here we try to solve the 
problem by constructing martingale [34] and bring Azu- 
mas inequality [35j into statistical fluctuation analysis. 
Detailed proof is shown in Appendix C. Here the quan¬ 
tum version of Azurnas inequality is presented as follows. 

Let X denote the observed value, which is obtained 
by detection events of m independent random samples, 
M n = i • X)T=i and Mq = E, here E is the value of 
1- ■ 5D"=i when n tends to infinity. Then for any £ > 0 
and £ > 0 , the following inequality holds 


Pr (| M n - Mo| > 0 = Pr (| ± ^ ~ E | > 0 

< 2e _2E ?=i x 'i , 

(13) 

with £ = yj ■=■ ln 7 . Applying Azuma’s inequality in solv¬ 
ing the statistical fluctuation problems based on depen- 
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dent samples, we get the relevant deviations denoted by 
£5 (S obs ) in Appendix El 


VI. SIMULATION 

To obtain the impact on secure key rate of after-pulse 
contribtions, we take Wei’s protocol as an example to see 
the comparisons with and without after-pulse effect. 



FIG. 1. (Color online) Plot of secure key rate R versus total 
transmission loss t with and without after-pulse contributions. 
The blue solid line denotes the secure key rate under after¬ 
pulse impact, and the blue dotted line stands for the secure 
key rate without after-pulse contributions. Here we take N = 
6 x 10 9 , u a = 5. 

Figure Q] shows that the secure key rate under after¬ 
pulse impact is 75% lower than the one without after¬ 
pulse contribution. As a result, the influence of after¬ 
pulse on secure key rate can’t be neglected. 

Take the after-pulse contribution into consideration, 
we obtain the secure key rates under different statistical 
fluctuation analysis methods in figure [2] 

And the corresponding deviation comparisons are pre¬ 
sented in figure [3] 

From the above analysis, we get a clear comparison of 
the estimation accuracy using different methods. And 
for different relations of the observed samples, we should 
choose the proper way to deal with statistical fluctuation 
problems so as to get a precise result. 


VII. CONCLUSION 

In this paper, we give a detailed analysis in dealing 
with statistical fluctuation issues, based on independent 
and dependent random samples. We set up respective 
mathematical models for each estimation method and 
give a fair comparison standard, together with the cor¬ 
responding deviation comparisons. For the first time we 
put forward a stringent way in handling the statistical 



FIG. 2. (Color online) Plot of secure key rate R versus to¬ 
tal transmission loss t under different estimation methods. 
The red lines denote the respective secure key rates with 
N = 5 x 10 10 , and the blue lines refer to the ones with 
N = 10 12 . The abbreviations respectively stands for Law 
of large numbers (LLN), standard error analysis (SEA), Ho- 
effding’s inequality (HI), Chernoff bound (CB) and Azuma’s 
inequality (AI). Taking N = N to tai £ {5 x 10 lo ,10 12 }, e = 
£1 = £2 = £3 = 10~ 10 , and u a = 6.4. 



FIG. 3. (Color online) Plot of deviations comparisons versus 
total transmission loss t under various estimation methods. 
The blue lines from top to bottom respectively stands for 
the deviations of Law of large numbers, Azuma’s inequal¬ 
ity, denotes the deviations when using law of large num¬ 
bers, Azuma’s inequality, Hoeffding’s inequality, Hoeffding’s 
inequality , Chernoff bound and standard error analysis with 
N = 10 12 . 


fluctuations based on dependent samples. From the re¬ 
sult of numerical simulation, we conclude that the after¬ 
pulse impact cant be neglected since it greatly lowers the 
secure key rate. And if the number of the totthe estima¬ 
tion accuracy of standard error analysis is the best one 
in the four methods. However, the large deviation the¬ 
ory possesses more generality than standard error anal¬ 
ysis. Although the secure key rate using Azumas in¬ 
equality is a little lower than Hoeffdings inequality, its 
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a strict tool dealing with fluctuation problems when tak¬ 
ing after-pulse contributions into consideration, and can 
be extended in other QKD protocols. 

We emphasize that Azurnas inequality is a stringent 
way in dealing with statistical fluctuations based on de¬ 
pendent samples, but not an efficient one. So there is still 
more room for accuracy improvement in finding other 
proper ways. 


Appendix A: Respective deviations 


Deviations under standard error analysis are repre¬ 
sented as follows. 


6 (Q„) 

Ci (Qv) 
Ci (bo') 
Ci (Qo) 


y/NpPzQfj, 

Qv'U’OC 

y/NjjPzQv 


Ug 

VNoQo 


(Al) 


follows. 

Ql = Qu (! + C 4 (Qij.)) = Q» + \J ^7 ln ^ 

Qv = Qv{i-U(Qv)) _ 

= Q v - min In A, ^2s_ In 

y 0 l = n(i-c 4 ( v)) _ 

= Y ° ~ min {v/^ In b \/^r In Tb } 

Qo = Qo (i - Ci (Qo)) 

= Qo- mm l^/^lii A, \J In } 

(A4) 

Deviations under Azuma’s inequality are represented as 
follows. 


Q^ — Qn (1 + Cs {Qij,)) — Q/j + y -jy-~ In 7 

Qv = Qv (1 + C 5 (Qv)) = Qv — \ ITfZ In — , 

y A r~p z s (A5) 

Y 0 L =Y'( 1 + Cs (bo')) = bo' - \j^n~ £ 

Qo = Qo (1 + Cs (Qo)) =Qo- \j-jy- In - 


where Q M , Q v , boand Qo are the observed values instead 
of probabilities. 

Deviations under Law of large numbers are represented 
as follows. 


Appendix B: Chernoff bound 

When it comes to Chernoff bound theory, there are 
three conditions satisfied or not corresponding to six dif¬ 
ferent results. 


Ql — Qic (1 + C2 (Qm)) — Q/j + 


2 [ln(b) + 21 n(jy z + l)] i. ( 2 ^ 1 ) 1/roat < e (3 / 4 ^) 
N^p z ii. (e^) 1/maL <e^ 

Qv = Qv (1 - c 2 (Qv)) = Qv - + 4 iii.^ 3 1 ) 17 - < ~ 

V _ *vPz _ Let f(x.v) = ./^hU. 


Y 0 l = Yo'( 1 - C2 (Yo 1 )) = Yo' - 


Let f (x,y) = W-^ln-, then we have the following 
2[ln(I)+21n(iVo + l)] results . ' V ™ " 


No 


(1) If condition (i) and (ii) are fulfilled, then A = 


pL r a C ,p p / 2 [ln(i) + 21 n(A 0 + l)] /(*,ef/l6) and A = / (x, e 3/2 ). 

Qo - Qo (1 - C 2 (Qo)) - Qo - y-]Vo- ( 2 ) If condition (i) and (iii) are 


(A2) 

Deviations under Hoeffding’s inequality are repre¬ 
sented as follows. 


fulfilled, while (ii) not, 
then A = f (x, £ 2 / 16 ) and A = / (x, £§). 

(3) If condition (i) is fulfilled and (iii) not, then A = 
f(x,4/ 16) and A = f(m/A,e). 

(4) If condition (ii) is fulfilled and (i) not, then A = 


Qfv — Qfi (1 + C3 (Q/j)) ~ Qfj + 

Qv = Qv (i - Cs (Qv)) =Qv-\ 


—-— In — 

2N^p z e 

1 ■ In - 


/ (m/ 4, e) and A = / (x,i 


.3/2^ 


Yo = Yo (1 - Cs (bo')) = Y o ~ \j 
Qo = Qo (1 — C3 (Qo)) = Qo — 


2 N*p x £ (A 3 ) 


— Ini 

2Ap £ 

1 , 1 

2 V 0 11 £ 


Deviations under Chernoff bound are represented as 


(5) If condition (i) and (ii) are not fulfilled, while (iii) 
is, then A = / (m/4, £) and A = f (x, £§). 

( 6 ) If condition (i) , (ii) and (iii) are not fulfilled then 
A = A = / (m/4, e). 

The other deviation comparison results with Hoeffd¬ 
ing’s inequality are stated here, also taking Q M as an 
example and £1 = £2 = £3 = 10 -1 °. 

(a) Deviation comparison between yj ln If, 

\/ Sr ln ^ and y/ 2i^: ln i- under result ( 2 ) - 
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The sufficient condition that Chernoff bound achieves 
a better result than Hoeffding’s inequality is 


'-^ln“< 

NfiPz £ 4 


I 2 


In ■ 


1 


.2 < 


'—— ln- 

2 N^pz e 

CZ^I 

2I+p z £■ 


(Bl) 


This condition can be equivalently written as 

<+ < 0.06. (B2) 

(b) Deviation comparison between ^ In and 

^ 2iv^r ln b under result ( 3 )- 

Obviously, the condition is the same as (a). 

(c) Deviation comparison between ln -j— and 

\J 2W~p^ ln e ’ UIlder reSult ( 4 )- 


Here the sufficient condition is 


/ 2Q P 

NfzPz 


72 


—-—ln — 

ZN^pz £ 


And we have 


Qn < 0.162. 


(B3) 


(B4) 


(d) Deviation comparison between yj ln and 

\J 2lV^T ln e"> Under reSult ( 5 )‘ 


Here the sufficient condition is 





(B5) 


And we have 


< 0 . 122 . 


(B6) 


Next we construct the martingale through analyz¬ 
ing the mutual relationships between dependent sam¬ 
ples. If the observed value X are obtained by detec¬ 
tion events of in dependent random samples denoted by 
X \, X 2 , ■ ■ ■ ,X m , detected with the value 1 and 0 other¬ 
wise, satisfying X = X Xj,. Considering after-pulse 
contributions we have the following probability relations 


Pr s| X m — t, X n —\ — x n —i 5 , Aq — 3 ^ 0 } 

= Pr {X m+1 = s\X m = t} 

£ Pr {X m+1 = s} 

(C2) 

(D2) shows that the relationship of the samples sat¬ 
isfies a Markov chain. And take signal pulse as 
an example, we have P {X m+ i = 1} = Qp m and 

P {-^m +1 — 1 |Am 7er( 1 Qp m and Ep m 

are the measured values of overall signal state gain 
with and without considering after-pulse contributions, 
when the number of observed samples is m, satisfying 
•D Mm = Qvm^+Pap)- Let S n = Yh= i X i = then 
= — — — • And we have the following relationship. 

E(S n+ 1 1 S n = k) 

= k ■ P(S n+ 1 = k\S n = k) 

+ (fc + 1) • P(S n+ i = k + 11(S'rj = k) 

= k ■ P{X m +i = 0 \X m = t} /po\ 

+{k + 1) • P{X m+1 = l\X m = t} [Kj6) 

= k ■ (1 — D/j, n ) + (k + 1) • Dfz n 

= k + D ltn 
= S n + D lln 

Let M n = ^,Mq = E , here E is the value of L • 
y^7-i Xi when n tends to infinity. We have 


E(M n+1 \S ir -- ,S n ) 


E{M n+1 \S n ) 

Sn) 

4 E(S n+1 \S n ) 
1 (S n + D^) 


n+1 

S n +S. 


V™ 


n+1 

1 = Mr, 


(C4) 


(e) Chernoff bound changes its form into Hoeffding’s 
inequality under result (6). 


From the above analysis we can conclude that M n is a 
martingale which can be used in Azuma’s inequality. 


Appendix C: Detailed proof 

First let’s look at the definition of Azuma’s inequality. 

Definition:Let {Y 1 ,Y 2 , ■ ■ ■ , Y n } be a martingale, and 
suppose that there exists a sequence {X 1 ,X 2 ,- ■ ■ , W n } 
of real numbers such that Pr (| Y n — W-i| < K n ) = 1 for 
all n. Then 

Pr(|y„ - Fol > 0 < 2e ( -^ 2 /^ x '),e > 0. (Cl) 
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